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TEHD lubrication of a cylindrical and elliptical 
bearing two lobes "parametric study" 
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Abstract — In this article, the thermoelastohydrodynamic 
study for the analysis of cylindrical and elliptical bearing 
(two-lobe) with a Newtonian lubricant was presented. The 
thermoelastic deformations of the solid parts are taken into 
account. To solve the Reynold's equation generalized, equation 
of energy and the displacement field, respectively, using two 
numerical techniques Computational Fluid Dynamic "CFD" 
and Fluid Structure Interaction "FSI". This article will aim to 
compare a cylindrical and elliptical bearing with two lobes. This 
comparison will focus on the various operating parameters of a 
bearing (pressure, temperature, displacement and stress) for 
different values of speed and eccentricity. 

Index Terms — Thermoelastohydrodynamic, Cylindrical 
Elliptical journal bearing, Two-lobe. 


2010 [3], conducted a study "THD" parameter for the 
analysis of behaviors of an elliptical bearing lubricant for 
three different category. The finite difference method was 
adopted for the numerical solution of the generalized 
Reynolds equation and the energy in the lubricant film. Vikas 
Phalle, Satish C. Sharma and SC Jain [4] (2011), in this work 
the authors present an analytical study on the influence of 
wear on the performance of an offset two -lobe membrane 
four a hydrodynamic bearing pockets. The generalized 
Reynolds equation governing the flow of lubricant in the 
clearance space of a bearing system was solved using finite 
element method with Newton-Raphson. The numerically 
simulated results clearly indicate that the performance of the 
bearing is largely affected by wear. 


I. INTRODUCTION 

In the literature, we examined the simple case of cylindrical 
bearing, in practice, other bearing types are often used as the 
cylindrical bearing is stable only for large eccentricities, that 
is to say high loads. Another type of fixed geometry bearing 
and has good dynamic characteristics is the lobed bearing. Its 
configuration is similar to that of the tilt pad bearing, except 
that the pads (lobes) are fixed in a given position beforehand. 
If the tip is angled so that the minimum thickness of the film 
is located at the outlet of this lobe, while the converging zone 
and the lifts are increased. The basic principle of non-circular 
bearing is to make changes in the symmetry of a plain 
bearing, which encourages the characteristics of self-excited 
vibration, and overcome the adverse effects of the 
configuration of circular rotation. To combat the dynamic 
instability of rotating at high speed in their bearing trees is 
often chosen non-circular plain bearings. Among these, the 
elliptical bearing has a certain predilection. 2007 K 
Prabhakaran and all [2] showed static and dynamic analysis 
of behavior "EHD" an elliptical bearing (two lobes) with the 
use of two types of Newtonian lubricant and microfleece. The 
finite element method is used to solve the modified Reynolds 
equation and the equations of three-dimensional elasticity 
directing the displacement field in the bearing. The static and 
dynamic behavior characteristics are calculated for a wide 
range of deformation coefficient that takes into account the 
flexibility of the bearing. The calculated results show that 
increasing the concentration of additive volume and mass 
transfer additive produced significant on the performance 
characteristics of the bearing change. Amit Chauhan et al in 
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II. EQUATIONS OF THE PROBLEM TEHD 
A. Lubricant film thickness 

For an elliptical bearing the lubricant film thickness is written 
as follows [3]: 



Figurel: Schematic elliptical journal bearing [3] 
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B. Reynold's equation 

The study of the lubrication of a bearing requires the 
application of various assumptions about mechanical viscous 
thin films. Assuming the fixed bearing and the shaft of a 
lively angular velocity co, writing this equation in a Cartesian 
coordinate system leads to the following formula [5]: 

d , h 3 5P „ 2 3 , h 3 dP s UR dh 

— ( ) + /?—( ) = 

d6 12// dO dz 12// dz 2 dO 

V_ h dU_ dh 

+ R 2 + 2R 2 do dt 

The boundary conditions on the pressure are quite simple: it 
is assumed that the reference pressure is atmospheric 
pressure. The pressure in the inlet sections and outlet of the 
film is considered equal to the supply pressure and the 
bearing edges of equal to atmospheric pressure. 

C. Energy equation 


The energy equation in mechanics of viscous thin films is a 
particular form of shape differs from an author to another 
depending on changes made to it. The energy equation 
written for an incompressible Newtonian fluid and the 
simplified calculation domain as follows. 
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Viscous dissipation 

(u, v, w): Component of velocity, 

[i\ Dynamic viscosity of lubricant; 

C p : Specific heat of lubricant, 

K h \ Conduction coefficient. 


characteristics of the bearing and the operating conditions 
and the parameters of the lubricant specified in Table 1. 


Bearing length 

Leu = 75mm 

Bearing length 

L cy i = 75mm 

Length/Diameter 

L C y/D cy i=0.57 

Outer radius of the bearing 

R ext = 72 mm 

Inside radius of the bearing 

R int = 65mm 

shaft radius 

R a = 45mm 

Young 's modulus 

n 

tNj 

+ 

2? 

Poisson coefficient 

v c = 0.33 

Thermal conductivity 
of bushing 

K hc =50W/m.K 

Thermal conductivity of shaft 

K ha =75W/m.K 

Lubricant viscosity 

Ho = 0.02068Pa.s 

Density of lubricant 

p=869. 53 Kg. in 

Specific heat of lubricant 

C p =2000J.Kg.K 

Thermal conductivity 
of lubricant 

K h =0.13W/m.k 

Supply pressure 

P o=l MPa 

Supply temperature 

T 0 =300K 

Rotation speed 

1 00 0<co< 4000rpm 

Eccentricity 

20pm< e <8 Opm 


TableI: The geometrical characteristics of the bearing 


(cylindrical and elliptical), lubricant and operating conditions 


Relative difference 

Convergence criterion 

A P/P 

I0 e 6 

A77T 


A h/h 


A u/u, Av/v, A w/w 

Kf 1 


D. Displacement equation 


TABLE2: Convergence criteria 


The equations governing the displacement field are given as: 

[M,]{(7}+[^]{£7} = [F,]+[R]{P} (4) 


M, 0 
pRT M f 



K S ~R S 

0 K f 




(5) 


Where, [ M s ] is the structural mass matrix; [Mf is the fluid 
mass matrix; [F s ] and [F f ] is the structural and fluid force 
matrix; [R] is a coupling matrix that represents the effective 
surface area associated with each node in fluid structure 
interface. The interaction of the fluid and the structure at the 
mesh interface causes the pressure to exert a force applied to 
the structure and the structural motions produce an effective 
"fluid load". 


III. Results and discussions 

The example given relates to a landing with two elliptical 
lobes, the pad in Bronze and Steel shaft. The geometrical 


A. Results CFD and FSI 



Figure 1: Mesh film lubricant (11910 nodes) elliptical 
bearing. 
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Figure2: Mesh of bearing (5200 nodes) elliptical bearing. 


Figure3 respectively each have the three-dimensional field 
distribution of pressure and temperature in the lubricating 
film for a rotation speed of tree 2000rpm and an initial 
eccentricity 20pm. It is observed, that the fields of 
temperature and pressure is commonly divided into the 
circumferential direction, and they reach its maximum axial. 
Thus respectively order of 18.2MPa and 330K. More 
pressure and temperature attained its maximum in the lobes 
of the bearing (0=-9O° and 0=90°). 
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Figure3: The 3D field distribution of temperature and 
pressure in the lubricant film 
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Figure4: The 3D field distribution of displacement and stress 
of the bearing 


Figure4 have the three-dimensional distribution of 
displacement fields respectively Usum and Ssum in the bush. 
Under the effect of the field of pressure and temperature, 
deforms the pad, to a greater extent in the median plane of 
which the thickness of the lubricant film is minimal. It may 
be noted that the stress peaks in the downstream of the area 
where the lubricant film thickness is minimal, which 
corresponds to the area of displacement Usum Maximum 
pad, especially at the ends Z=+L/2. The displacement field 
peaked in the lobes, order Umax= 68,4pm. The stress reaches 
its maximum value Smax= 87.1 MPa in the ends Z=+L/2. 
This is due to made that, at these ends; the installation (a 
boundary condition) generates significant constraints. 

IV. Parametric study 

In this part we are interested in studying the influence of 
the rotational speed on the field of pressure, temperature 
displacement and stress intensity in the middle (z=0 and 
O°<0<36O°). 

A. The influence of the speed on performance TEHD 


Figure5 expose the influence of the rotational speed on the 
field of pressure and temperature as a function of the angular 
coordinate in the median plane of the bearing. An increase in 
pressure around 5.29% of CDi=1000rpm and co 4 =4000rpm and 
an increase of order 8.22% for temperature is observed. 
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Figure5 : Circumferential distribution of fields P and T for 
different values of co 
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Figureb: Circumferential distribution of fields U and S for 
different values of co 

Figureb shows the influence of the rotation speed of the 
circumferential distribution of the displacement field U and 
the stress S in the bearing inner surface. For a speed of 
lOOOrpm was Umax=34.3pm and Smax=20.5 MPa and 
4000rpm found Umax=8bpm and Smax=81MPa therefore 
the displacement field and stress acceded respectively an 
increase of around 60% and b3%. 


B. The influence of eccentricity on performance TEHD 

In Figure7 we see that there is an asymmetry in the field 
distribution of pressure and temperature. For an eccentricity 
20pm pressure and temperature reached their maximum 
values Pmax=25.9MPa Tmax=308K and an eccentricity 
80pm pressure and temperature to Pmax=b8.1MPa 
Tmax=341K by finding the two variables respectively 
underwent an elevation of about 38% and 10.41%. 




Circonferentielle direction 


Figure7 : Circumferential distribution of fields P and T for 
different values of eccentricity 
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Figure8: Circumferential distribution of fields U and S for 
different values of eccentricity 

The influence of eccentricity on the displacement field U and 
the constraint S is presented respectively in Figure8. 

The displacement field has reached a maximum value of 
9b. 20pm for eccentricity 20pm and eccentricity 80pm 
reached a maximum value of 744pm, an increase was 
observed significant order 87.07%. Regarding stress, there is 
an increase of the maximum order value 70.50% between the 
eccentricity 20pm Smax=3b.2MPa and 80pm 
Smax=llbMPa. 



AP(MPa) 

AT(K) 

AU(MPa) 

AS (MPa) 

co 1= lOOOrpm; 
e3=b0pm 

2.55 

17,69 

111 

235 

col=2000rpm; 

e3=b0pm 

12.56 

27,99 

131 

274 

col=3000rpm; 

e3=b0pm 

27.6 

44.98 

149 

326 

col=4000rpm; 

e3=b0pm 

48,2 

56.98 

165 

351 


Table3: The influence of speed “co” on strained AP(MPa), 
AT(K), AU(MPa) and AS(MPa) a cylindrical bearing. 
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AP(MPa) 

AT(K) 

AU(MPa) 

AS(MPa) 

ol=2000rpm; 

e3=20|jm 

8.34 

18.12 

46.2 

94.77 

ol=2000rpm; 

e3=40)im 

9.22 

22 

54.7 

100.32 

ol=2000rpm; 

e3=60|jm 

12.65 

27.99 

131 

273.83 

col=2000rpm; 

e3=80jjm 

17.4 

34 

154 

297.67 


Table4: The influence of eccentricity “e” on strained 
AP(MPa), AT(K), AU(MPa) and AS(MPa) a cylindrical 
bearing. 



AP(MPa) 

AT(K) 

AU(MPa) 

AS(MPa) 

ol=1000rpm; 

e3=20|jm 

6.37 

15.57 

34.3 

19.4 

ol=2000rpm; 

e3=20|am 

27.68 

25.47 

38.5 

28.1 

ol=3000rpm; 

e3=20|jm 

54.17 

41.51 

68.6 

51.9 

O)l=4000rpm; 

e3=20|jm 

77.10 

46.70 

86.0 

78.4 


Table5: The influence of speed on strained AP(MPa), 
AT(K), AU(MPa) and AS(MPa) an elliptical bearing. 



AP(MPa) 

AT(K) 

AU(MPa) 

AS(MPa) 

col=2000rpm; 

e3=20|jm 

27.68 

25.47 

38.5 

28.1 

ol=2000rpm; 

e3=40|am 

54.71 

32.86 

45.1 

50.26 

col=2000rpm; 

e3=60|am 

79.04 

39.78 

60.9 

96.7 

col=2000rpm; 

e3=680|jm 

92.19 

50.7 

74.4 

115.29 


Table6: The influence of eccentricity “e” on strained 
AP(MPa), AT(K), AU(MPa) and AS(MPa) an elliptical 
bearing. 


Presenting in Figures 9, 10, 1 1 and 12 curves that allow us to 
provide an overall summary of the influence of the rotational 
speed and eccentricity variations on extended AP(MPa), 
AT(K), AU(MPa) and AS(MPa) for two geometries treated, 
cylindrical and elliptical. See Tables 3, 4, 5 and 6. 



Figure9: The influence of the rotational speed and 
eccentricity of AP(MPa) for both elliptical and cylindrical 
geometries. 



Figure 10: The influence of the rotational speed and 
eccentricity of AT(MPa) for both elliptical and cylindrical 
geometries. 



Figure 1 1: The influence of the rotational speed and 
eccentricity of AU(MPa) for both elliptical and cylindrical 
geometries. 



Figure 12: The influence of the rotational speed and 
eccentricity of AS (MPa) for both elliptical and cylindrical 
geometries. 
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Whatever the operating conditions and the geometrical 
characteristics of the bearing, it can be concluded that the 
hydrodynamic performance trends obtained through the two 
geometries are the same. 

Figure9 shows that the pressure values at the elliptical 
bearing are greater than those of the cylindrical bearing, for 
different values of and the eccentricity of the rotation speed. 
While the temperature values for the elliptical bearing are 
less than those of the cylindrical bearing figure 10, for 
different values of eccentricity and the rotation speed. 
Despite the significant values of the pressure field for the 
elliptical bearing compared to the cylindrical bearing, we see 
that the displacement values at the elliptical bearing for 
eccentricity 20pm and with different values of speed are 
lower. This shows the advantage of the geometric shape of 
the elliptical bearing when the rotational speed varies 
figurell. Unlike the influence of rotational speed on the 
displacement, that shows the influences of the eccentricity on 
the maximum displacement values are greater for the 
elliptical bearing. 
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V. Conclusion 

This paper has studied the nonlinear analysis of lubrication 
"TEHD" a cylindrical and elliptical bearing (two lobes). 
Using one hand, the numerical simulation Computational 
Fluid Dynamic “CFD” which allowed us to determine the 
pressure field and temperature in the lubricant film. On the 
other hand we used another Fluid Structure Interaction 
numerical simulation "FSI" to calculate the displacement 
field and the stress in the bearing. It should be noted, 
however, that the conclusions we give in this article are valid 
only for the specific cases that we studied, and they are not 
independent of the characteristics of the bearing and 
lubricant. Our work must be completed to show the influence 
of material properties and lubricant as well as the dimensions 
of the bearing elements in determining its performance. 
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